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Abstract. The Gross-Pitaevskii equation with a local cubic nonlincarity that describes 
a many-dimensional system in an external field is considered in the framework of the com- 
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direction of the surface normal. The general constructions arc illustrated by examples. 
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1 Introduction 

The Gross-Pitaevskii equation (GPE), derived independently by Gross and Pitaevskii [2], 
arises in various models of nonlinear physical phenomena. This is a Schrodinger-type equation 
with an external field potential U(x,t) and a local cubic nonlinearity: 

-ihd t + ^- + U{r) + x\V{x,t,h)\ 2 ) V(x,t,h) = 0. (f) 

Here x = (xj) £ W 1 ; i, j, k, I = 1, . . . , n; t G R 1 ; dt = d/dt; p = —iftV; V is a gradient operator 
in x\ x and h are real parameters. The complex function *$>(x,t,h) determines the system 
state, |^| 2 = t, h)\ 2 = *$>(x,t,h)fy*(x,t,h); the function $>*(x,t,h) is complex conjugate to 
$?(x,t, h). 

For U = 0, xr < 0, and n = 1, equation is known to be integrable by the Inverse Scattering 
Transform (1ST) method [31 15 an d to have exact soliton solutions. The solutions are written in 
terms of space localized functions which conserve their form during the evolution. Solitons have 
numerous physical applications; for example, they serve to describe the propagation of optical 
pulse in nonlinear media [SI El • 

The solutions of equation (^Q) that are not completely localized are also of interest, in parti- 
cular, in the theory of nonlinear deep-water waves where the two-dimensional (n = 2) solutions 
of plane wave solitonlike type are localized only along the wave propagation direction [7] . 



2 



A. Borisov, A. Shapovalov and A. Trifonov 



The Gross-Pitaevskii equation Q in physical dimensions (n = 2,3) is used in the mean- 
field quantum theory of Bose-Einstein condensate (BEC) formed by ultracold bosonic coherent 
atomic ensembles (Hj. As a rule, the function U(x,t) is the potential of an external field of 
a magnetic trap and laser radiation. The wave function ^>(x,t,h) corresponds to a condensate 
state. The local nonlinearity term x\^S>(x, t, h)\ 2 arises from an assumption about the delta- 
shape interatomic potential. For multidimensional cases, no method of exact integration of the 
GPE is known, to our knowledge, and equation is usually studied for some special cases by 
computer simulation. However, this approach has natural restrictions and also has nontrivial 
specific features as long as localized solutions of equation (^Q) are known to be unstable and to 
collapse when n > 1, U = 0, and h < Therefore, to develop methods for constructing 
analytical solutions of the GPE (pQ) for many-dimensional cases is critical. 

In this paper we construct a class of asymptotic solutions for the (1 + ?i)-dimensional GPE 
with a focusing local cubic nonlinearity 

L(V)(x) = \-ihd t +H(f>,x,t) - g 2 \^(x,t,h)\ 2 ] V(x,t,h) = 0. (2) 

Here g is a real nonlinearity parameter and h is an asymptotic parameter, h — ► 0. The pseudo- 
differential operator Tt(p,x,t) is Weyl-ordered JU] and its symbol Tt(p,x,t) is quadratic in 

l ~ l 



where (a,b) = ^ aibi; the real functions — the (re x n)-matrix H. pp (t), the vector 7i,(x,t), and 
i=i 

the scalar Tlo(x, t) — smoothly depend on their arguments. These functions model the external 
fields imposed on the system described by equation 

A formalism of semiclassical asymptotics was developed for the GPE with a nonlocal nonlin- 
earity in [111112] . This equation was named the Hartree type equation. The semiclassical asymp- 
totics of the nonstationary Hartree type equation were also studied ^3 El El El El El El HO] • 
This formalism is based on the WKB-Maslov complex germ method, the idea of which was first 
put forward in |5J and then comprehensively treated in |22L 1231 124j . 

Further developments of the Maslov method |25l I2(ij (see also [^7j and references therein) 
allow one to construct the localized solutions of the Cauchy problem in the class of trajectory - 
concentrated functions. 

Here we apply the ideas of the complex germ method and the approach used in ^2 Ej to 
construct the analytic solutions asymptotic in a small parameter H, h — > 0, for equation P|). 
The asymptotic solutions are constructed in the class of functions localized in the vicinity of 
an unclosed surface associated with the phase curve that describes the evolution of the surface 
vertex. 

The functions of this class have the one-soliton form of the one-dimensional nonlinear Schro- 
dinger equation along the surface normal direction. A semiclassical linearization of equation © 
is performed accurate to 0{ti i / 2 ), h — > 0. The asymptotic solutions constructed are illustrated 
by examples. 

2 The class of paraboloid-concentrated solitonlike functions 

Let us describe a class of functions ^(x,t,h) in which the asymptotics for the GPE © are 
constructed. Following [111 I12j . consider equation ^ for the (1 + l)-dimensional case with no 
external field (U = 0) 

(ihdt + hi 2 d 2 x + g 2 \^(x,t,h)A *(x ) t,h) = 0, (4) 
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where iGl 1 and d x = d/dx. Equation © has a one-soliton solution [3J|1] 

(5) 



2r\ 1 

^(x, t,H) = —. — ; . 7-z exp 

v ' g ch[2??(x-xo -2£t)/h] 



l -(2tx-2{e-r, 2 )t + w) 



where £, rj, xq, and ip® are the real parameters. 

Earlier |28| we obtained a semiclassical solution of equation (J2J) with the leading term 

t, K) = g-^iVa^ppi^Va) exp Us(x, t, ft) J ch" 1 t, &)) • (6) 

It was assumed that the complex function Q(x,t,h) = S(x,t,h) + ia(x,t,H) regularly depends 
on h and its high-order part, i.e., the function $l(x,t, 0) = fl,(°\x,t) = S^°\x,t) +i&( '(x,t), is 
a complex solution of the Hamilton-Jacobi equation 

d t n {0) +n{vn {0 \x,t) =o. (7) 

The functions © are localized in the neighborhood of a i-parameter family of surfaces 

r* = {x : a(x,t,0) = a (0) (x,t) = 0}. (8) 

Assume that the rank of Hesse matrix of the function a^°\x, t) is n at any point x G T*. 

The construction of global solutions of equation (JJJ) is beyond the scope of this work. Here 
we consider the evolution of the functions localized in the neighborhood Ug of the surface 
point xo © at the initial time t = 0. 

Denote by 

7 = {x : x G M n ,x = X(t),t ^ 0} 

a smooth curve passing through the point xq (= X(0)) in the coordinate space x G R n such that 
-X'(i) G r . This curve plays the role of the "classical trajectory" corresponding to the solution 
of the quantum equation ©. 

In the neighborhood Uxt t \ C T* of the point X(t) G 7, the hypersurface T* can be approxi- 
mated by a simpler surface, for example, by the tangent plane. It is well known that a tangent 
plane is determined uniquely by its normal and the point of contact. Let 7?(t) be a normal vector 
to the surface T* at a point x = X(t). 

To construct asymptotic solutions of equation © one needs the complex germ |22| . i.e. the 
n-dimensional complex space associated with the equation under consideration and functions 
defined in the neighborhood of the hypersurface (JSJ). 

It is more precise to use a second order surface 

(7?(t), Ax) + -(Ax,Q 2 (i)A£) = 0, Ax = x-X(t), (9) 

instead of the tangent plane. The surface @ with a proper choice of the matrix Q%{t) is 
a quadratic approximation of (JSJ) in the neighborhood C^jf ^ • 

Introduce the class N\ of complex functions with the generic element 

*(x, t, K) = $(0, x, t, h) = (T 1 ^{^(t),Hpp(t)Tr(t)) exp Us(x, t, ft) J 

x ch -1 (0){i + hch(e)u(e,x,t,n) +ihch(e)v(e,x,t,h)}. (10) 

The functions (|1U|) are localized in the neighborhood ^x(t) of the point x = X(t). This point is 
the "vertex" of the surface ©• In equation ()l(Jj) . = a(x,t,H)/H is a "fast" variable; the real 
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Figure 1. 



functions u 
in 9, and 



, x, t, h) and v = v(6, x, t, h) are regular in their arguments and are bounded 



a(x,t,H) = (7f(t),Ax) + -(Ax,Q 2 {t)Ax) + ha 1 (x,t,h), 
S{x,t,h) = S°(t) + (P{t),Ax) + i(Ax,Qi(t) As) + hS 1 (x,t,h). 



(11) 
(12) 



In equations (|1 If) and (|12|) . the real function S°(t), the real n-dimensional vector functions 
Tr(t), P(t), and X(t), the real (n x n)-matrix functions Qi(t) and Q2(t), the real functions 
cr 1 (s, t, h) = cj 1 (^, i, h) and S^^r, i, ft) = S 1 ^, t,h), £ = Axj\ffi are functional parameters of the 
class JVJ[; Af = x- X(t). 

Note that in constructing the class JV^, a specific basis, "orthogonal" in some way to the 
vector 7r(t), is used 1 . As a result, the classical trajectory depends on the complex germ as well. 

The functions (jlOj) are not normalizable in the space L2(M. n ,dx) since they are localized on 
an unclosed surface, equation ©. Therefore, the surface r* (JHJ) (compact in some physical ap- 
plications) is substituted by its quadratic approximation Q in the neighborhood Ug^. As we 
are interested in finding functions localized in Ug, t y we come to the problem of small pertur- 
bations on the "background" of the hypersurface (jyj). To eliminate the "background", consider 
the function 



*(x,t,h) = ${9,x,t,h) - *(9 ,£,t,h) « g ~ l — ^(K,H pp (t)K) chr\9) 

x exp (?-S(x,t,h)j {1 + hch(9)u(9,x,t,h) + ihch(9)v(9,x,t,h)} 



01, (13) 



-(7f(t),A5> + — (Ax,Q 2 (t)Ax), 9 1 = a 1 (i,t,h). 
n 2n 



The functions *$>(x,t,h), equation are normalizable in L2(W l ,dx). Estimating the so- 
lutions of equation P|). we use the norm of the functions *$>(x,t,h). Taking into account this 
normalization condition, let us assume that the functions a 1 (^,t, h) and S' 1 (^,t, K) belong to 
the Schwartz space 5(R^ _1 , dfi), where dfi is the measure on the hypersurface, determined by 
equation ©. 

The class J\fi of the form (|10|) with n = 1 includes the exact one-soliton solution © of one- 
dimensional nonlinear Schrodinger equation @ with a special choice of the functional parameters 

1 Such a situation is typical of the complex germ method 1221 1241 where the solution of spectral problem is 
considered for the functions localized on incomplete Lagrange manifolds with the use of a germ basis which is 
"tangent" to the manifold. 
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of the class. If the rank of matrix Q2(t) is n, then in one-dimensional case equation (110(1 can not 
be transformed to © since the argument 9 of the hyperbolic cosine contains the terms quadratic 
in x. The sufficient conditions for such a transformation are 



rank (Q 2 (t)) = n — 1, rank (Q 2 (*), = n, 



(14) 



where (Q2(t),Tr(t)) is the augmented matrix of order nx (n+1). Note that under conditions (ITU , 
the surface @ has the paraboloidal shape. 

Let us expand the operator L in equation (J2J in the neighborhood Ugr* under the conditions 



Ax k = d(h 1 / 2 ), Ap k 



6(ft 1/2 ), (?f(t),A2) = 6{h). 



(15) 



The leading term of the asymptotic in the class of functions (|1U|) under conditions (|15|) has the 
form ©. 

It is determined by the phase curve z = Z(t,H) = (P(t,H),X(t,H)), vector tt = n(t), and 
functions <j(x, t, K) and 5(x, t, K) of the form ((TT)) and (fT2|) . respectively. We call these functions 
the paraboloid- concentrated solitonlike functions. 



3 The linear associated Schrodinger equation 

The asymptotic leading term is obtained when the asymptotic solution of the equation (J2J) is 
constructed accurate to 0{h 2 ) 

Let us substitute (|1U|) in take into account the terms of order 0(H 2 ), separate the equa- 
tions with respect to the "fast" variable, and find the solutions of these equations that decrease 
as 9 — > ±oo. We than obtain for the function £1 = S + ia 



n t + H(vn,x, t) 

+ (n,H pp {t)7f) 



ih 

' y 

(it 



Sp[Wpp(t)^xa ; + U px {VVL, x, t)] 



{Tr,Hpp(t)n 



0. 



(16) 



With the substitution O = — ihln 

[-ihd t + H&x,t)]ip(x,t,h) = 0(h 3/2 ), 
where the operator 7i.(p, x, t) has the form of Q, 



i/j(x,t,H){Tr,7-lpp(t)jr} x l 2 . equation (fTB^l yields 



tp(x, t, h) = yj (jr,Hpp{t)Tr) 
ip{x, t, h) = {jr,H pp (t)ir) 



exp 
exp 



— u(x, t, h) 

n 



exp 



-S(x,t,h) 
n 



(p(x,t,h), 



(17) 



a(x, t, h) 



We call the equation 

[ - ihd t + H(fi, x, t)] ip(x, t, h) = 



(18) 



the linear associated Schrodinger equation corresponding to the GPE iffi). 

Therefore, the leading term of the asymptotic solution of the nonlinear equation @ is con- 
structed by using a solution of the linear equation (|18() with conditions (JBJ) and (|17|) . 

Note that for the function © to satisfy (|14|) the solution (|17|) must be chosen in a specific 
way. Denote the class of these functions by <S^. 
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Let us seek the solutions to equation (|TH|) that possess this property in the class of functions 



ip(x, t, h) = exp 



--(7f,Ax) 



(x, t, h). 



(19) 



Substitution of (JTHJ) in (JTHJ) yields 

[ - ihd t + i(if,Ax) - i(TT,x) + H{p{t) -\-m + Aft, x, t) ] cf>(x, t, H) = 0{ti i/2 ). (20) 

Let us expand the operator of equation (|20|) in a power series in Ax accurate to 0(h 3 / 2 ) in 
terms of the accuracy of ()15|) . We then have 

-ihd t -i[(n, (x--H p (t))) - (^ + n xp (t)TT),Ax) - (7r,n pp (t)Af)] 

+ H(t)- ^,H pp (t)n) + (Hg(t),Ax) + (H f (t), Ap) 



+ 



(Ax,H xx (t)Ax) + (Ax,H xp (t)Ap) + (Ap,H px (t)Ax) 



+ (Af>,n pp (t)Ap)]]^ + o{h 3 / 2 ) = [-ihd t + no(t)}(P + o{h 3 / 2 ) =o. 



(21) 



Therefore, for the solution ()19|) in the approximation under consideration, the linear associated 
Schrodinger equation (|T%|) takes the form 



[-ihdt+n (t)](f> = 0. 



(22) 



Here the operator TCo(t) (quadratic with respect to the operators Ax and Ap) is obtained 
from dZED . 

Thus, the construction of asymptotic leading term jf^) in the class of functions (|19|) for the 
nonlinear equation @ by solving equation (|22j) is complete. 

4 Solutions of the Gross— Pitaevskii equation 

The solutions of equation (|2*2*)) are well-known (see, e.g., EM ) • I n particular, the evolution 
operator U-]-t (t,s) is found in explicit form, whose action on the function ip(x,h) referred to 
a time t = s is given by the relation 



i){x, t, h) = U Ho (t, s)ip{x, h) = G Ho (x, y, t, s)ip{y, h)dy. 



(23) 



Here G-n (x,y,t, s) is the Green function of the linear associated Schrodinger equation (|2*2*|) . 
which is determined by the following conditions: 



{ - ihd t + Ho (t) } Gn (x, y , t , s) = 0, lim G Hq (x, y , t , s) = 5(x - y ) . 



(24) 



To construct a subclass of the class Mi (|1U|) of solutions of equation (J2J), let us find a particular 
solution in the form 



<f>(x,t) = iV^exp <{ - 



S°(t) + hS 1 ^) + iha l (t) + (P(t),Ax) + Uax, Q(t)Ax) 



(25) 



where the functions S°(t), S 1 ^), and cr 1 (t), the n-dimensional vectors P(t), X(t), and the 
complex (n x n)-matrix Q(t)=Qi(t) + iQ2(t) are to be determined; iV^ is a constant. 
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Substituting in (f2*H) and setting X = TC^, Q(t) = Qi(t) + iQ2(t), we obtain the deter- 
mining system of equations for P(t), 7r(i), X(t), Q(t) 

P = -Hg(P, X, t) + Q 2 n pp (t)7r, (26) 

x = n P iP,x,t), 

i? = -[n xp (P,X,t) + Qx(t)Hpp(t)]n, 

Q + H xx (t) + n xp (t)Q + QH px (t) + QH pp (t)Q = 0. (27) 
Consider the following system of equations 

Z = H px (P,X,t)Z + H pp (t)W, W = -H xx (P,X,t)Z -H xp (P,X,t)W. (28) 

Equations (|28|) are called a system in variations in vector form |27j . 

In general, it has n complex linear independent solutions, which can be written in the form 
of 2n-dimensional vector columns 

a j (t) = (W j (t),Z j (t)Y, j=T~H. (29) 

For equation (|22|) the vectors ()29|) set the symmetry operators which are linear in coordinates 
and momenta 

aj(t) = Nj [(£,-(*), A£) - (Wj(t), Ax)] , Nj = const, (30) 
since the conditions 

-ih^- + [H (t),a j (t)}=0 (31) 

are satisfied because of the validity of equation (|28|): [7^o(i)> = Ho(t)aj(t) — aj(t)7io(t) is 

the commutator of the linear operators. 

Let -B(i) and C(i) denote the (n x n)-matrices whose columns are constructed of the vectors 
of the solutions of system 1)2? 



B(t) = (Wi(i), W n (t)), C{t) = (Zi(t), . . . , Z n (t)). 

Then Q(t) = l?(t)C _1 (t) and = Q(t)Z(t), and the system in variations (|2*%)) is equivalent 

to (E3) (see EH). 

The normalization factors iV" are chosen as follows. The matrices B{t) and C(t) satisfy the 
condition 

Do(t) = hc + (t)B{t) - B + (t)C(t)} = A)(0) = A)- 
li 

Here C + is Hermitian adjoint to C. In addition, if the matrix Q(0) = -B(0)C _1 (0) is symmetric 
at the initial time t = 0, Q(0) = Q(0) T , then Q(t) = Q(t) T , t > 0. 
Let 

D = diag(d 1 ,d 3 ,...,rf n -i,0), dj>0, iVj 1 = (2M j )" 1 /2 ) N a n = l. 
Then the operators ()3U() satisfy the ordinary commutation relations 



[a,j(t),al(t)]=5 jk , j,k = l,n-l, 
and all the other commutators are equal to zero. 
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Let us denote the function ip(x,t, K) in lfTU|) as 



ip^\x,t, K) = iV exp 



n 



(7f(t),Ax> 



\o,t), 



(32) 



where the function (p(x,t,ft) of the form (|25|) is denoted as |0, t). The normalization factor iV"o 
is determined from the normalization condition for the function ^{x,t,ft) that is a solution of 
the nonlinear equation @. 

Let us substitute (J23) in $T^; then 



f detC(0) 
detC(t) 



exp 



■-(vf(t),Af) 



x exp 









IS 


[/'[ 

Jo . 



(P(t),X(t)) - H(t) + -(m,n pp (t)n(t)) 



dt 



+(P(t), Ax) + -(Ax, Q(t)Ax) 



Here (P(t), X (t)) is a solution of the first pair of equations (|2fi[). the (n x n)-matrix C(i) is 
composed of the columns Zj(t). Note in addition that (Ax, ImQ(t)Ax) > and 



aj(t)\0,t) = 0, j,k = l,n-l. 
We have for the function (x, t, ft) 



anitm {%,t,K) = 0, 



(33) 



(34) 



fi n (i) = exp 



-(vf(t),Ax) 



a n (i)exp 



--<7?(i),A5) 



a n (i) - i. 



Let l I / o(^) i> ft) denote the function ^>{x, t, h), Eq. (jHJ), corresponding to the function (x, t, h): 



1 (n,H pp (t)Tr) exp HSo(x, t, ft)) 



g 2 ch.(±a (x,t,ft)) 
In accordance with (|17[) . we have 



(35) 



S (x,t,h) 



(P(t),X(t)) - H(t) + ^(*),H pp (i)7f(t)) 



+ (P(t), Ax) + -(Ax, Qi{t)Ax) + Mm 



In 



'detC(O) 
det C(t) 



+ Mm N Q , (36) 



a (x,t,ft) = (n(t),Ax) + -{Ax, Q 2 {t)Ax) + ft In \f (n(t), H pp (t)n(t)) 



- ft Re 



In 



/detC(O) 
det C(t) 



ftReiY . 



(37) 



Equations (|35|) - (|37j) determine the leading term of the asymptotic solution of (j2j) that corre- 
sponds to the solution i/j^\x,t,h) of the form (|17|) for the linear associated Schrodinger equa- 
tion ((THJ). 
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Let us introduce a set of functions \v,t) with the help of "creation operators" cu(t), j 

(n(t),Ax) 



— acting on the function |0, t): 

1 



= N v exp 
= N v exp 



<7f(t),Af) 



1 



JJ |0,t) 



(38) 



Here v = 1/2, ■ ■ ■ , v n -i) £ is a mutiindex, a + (i) = (a^(i), . . . , a^_ 1 (t)), and H u (x, t, h) 

are many-dimensional Hermite polynomials. The normalization factors N v are determined from 
the normalization conditions similar to those for Nq in (|32l) . 
Using expressions similar to (|35j) - (|37|) . we obtain the function 



(vf, Hp P {t)T?} exp {{S u (x, t, h)) 



' g 2 ch (\a u (x,t,h)) 
S u (x, t, K) = S Q (x, t, h) + Mm [kiH v (x, t, h)] + Mm (N u - N ), 
a u {x, t, h) = (Jq(x, t, H) — HRe \hnH u {x, t, h)~\ — ft Re (N v — Nq) 



(39) 

(40) 
(41) 



for the function © corresponding to ^\x : t,K). 

Let us define class K, l h of solutions for the linear associated Schrodinger equation (fT%)) as 
a linear envelope of the functions ip^\x, t,h). The generic element of class /C| is given by 

tpic(x,t,h) = c u ^\x,t,h) = exp ( -i(vf(t), Ax) j (j} K {x,t,K) 
M=o V J 



kl=° 
The condition 

a n (t)ip)c(x,t,h) = a n (t)<f>ic(x,t,h) =0 



(42) 



is valid for the functions of the class K} h . This follows immediately from ()33|) . (|34(). and (|38() . 

Let 4>ic(x, s, h) be a function of the class /C| referred to a start time s, and t^ (i, s) is the 
evolution operator given by ()23j) and (|2"4*)) . 

Then the function 

1 



ip K (x,t,h) =exp I --(ir(t), Ax) \ U Uo (t, s)(f> K (x, s,ti) 

also belongs to the class /C^, since the function 

<f>K.(x,t,H) = U Ho (t,s)cf>ic(x,s,h) 
is represented as 

00 

(p K (x,S,h) = ^ C v\v,t). 
|i/|=0 



(43) 
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The above relation follows from the uniqueness of solution of the Cauchy problem for equa- 
tion (|22j). On the other hand, we have for the symmetry operator a n (t) of equation (|22j) the 
relation 

a n (t) = U Ho (t,s)a n (s)U^{t,s) 

following from (|23|) and 1)31 Jl. If the function <pic(x,s,h) belongs to the class /C| at the initial 
time s and satisfies the condition 

an(s)^ K (x,s,h) = 0, (44) 

then condition 1)42(1 is also satisfied for the function 1)43)1 at a time t > s, i.e., we have 

a n (t)4>/c(x,t,h) = 0. 

Therefore, the function (f)/c(x,t,h) belongs to the class K, l h . 

The evolution operator JJ-^ {t,s) of equation (|22|) induces the evolution operator U tr (t,s,-) 
for the nonlinear equation © in the class S^. 

Let a function tfj(x,h) be referred to an initial time s. According to (|17|) . the function 
^ s (x, H) corresponds to ip(x,h). Then the function ^(x,t,h), determined by (j^SJ), corresponds 
to ip(x,t,h). This correspondence can be written as a result of the action of the evolution 
operator Ut T (t, s, •) on the initial function ^ s (x, h), 

t>(x,t,h) = U tI (t,s,V s )(x,t,h). (45) 

Using the evolution operators U-H (t, s, •) and Ut (t, s, •) of the forms (f2"3*|) and (|43)) respectively, 
we can define the symmetry operators for the nonlinear equation (j2j) in the class of functions S^. 
To that end, let us take the function c/?(x, h) from the class K, s h which satisfies condition (|44j) at 
an initial time t = s. Let (p(x,t,h) be a function obtained from Eq. (|17[) and ^(x, t,h) be the 
solution of the nonlinear equation @ related to (p(x,t,h) according to (|17|) . 
Consider the operator A, A : K, s h — ► /C| such that 

[i,a n (s)] = 0, 

and a function (fA(x,h) = Aip(x,H) £ /C|. Then the function ^(x, t,h) related to ipA(x,h) 
by (|17|) can be treated as a result of the action of the symmetry operator A n \ of the nonlinear 
equation (J2J) 

^(x,t,n) = i nl ^(x,t,^). 

In conclusion, we note that for a nonlinear Schrodinger equation with a focusing nonlinearity, 
the many-dimensional solutions localized at the initial time are unstable. This leads to the 
phenomenon of collapse in the course of evolution. The semiclassical asymptotics (|39|) behave in 
a similar manner. They can be constructed for special external fields within finite time intervals 
where singularities typical of collapse appear. 

5 The three-dimensional anisotropic oscillator 

In case of a harmonic oscillator field, the linear operator 7i, equation Q, reads 
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where K = diag{/ci, k 2 , k 3 } with kj £ I 1 , j = 1, 3. The GPE © then takes the form 
P' 



^2 1 



-ihdt + ^- + ^(x, Kx) - g 2 \*(x, t, h)\ 2 
Zm Z 



$(x,t,h) = o. 



(47) 



To construct the asymptotic solutions (JBJ) for equation (|47|1 we solve the dynamic system (|26|) 
and (j28j) . which is reduced to 



-Kx H 7T, 

m 

1 



x = — p, 
m 



7T = 7T, 

m 



Z=—W, W = -KZ 
m 



(48) 
(49) 



in the case under consideration. Here Wj(t) and Zj(t), j = 1,3, are the linear independent solu- 
tions of the system in variations (|49|). which determine the matrices .B(i) = ( Wi(i), W 2 (t), Ws(t)) 
and C(t) = (Zi(t),Z 2 (t), Z 3 (t)); Q(t) = Qi(i) + iQ 2 (t) = B(t)C -1 (t), the matrices Qi(t) and 
Q 2 (*) being real. The solution of system (gSJ) and 103) is substituted in lf3"o]). (|37|) . (gDJ), and dHJ, 
and the functions S u (x, t, h), a u (x, t, K) are obtained. The function ^ u (x, t, K) is then determined 
by Pty. 

To write down the solution of systems (|48|) and (|49|) and the function ^(x, t, h), we introduce 
the notation fit = ^kj/m for kj > 0, QJ = y/—kj/m for fcj < 0, aj(t) = (Wj(t), Zj(t)) , and 
{a,j(t),ak{t)} denoting the skew-scalar product of the 6-vectors aj(t), ak(t), j,k = 1,3. 

Let us also introduce the functions 



Sn2 jn3 (t,X,p,Sli,n2,tt3,z(t)) = S°(t) + y^i^z(t) 



n 2 + - )n 2 + n 3 + - O 



7Tft , . . 

+ ^( n 2 + n 3 ) - ft 
o-na,^ (tto,*, ^ ^2, S2 3 , z(t)) = H — x 2 



t + (x,p), 



h In 



1 4 /?Tl 2 ^2^3 / 1 



a; 



"2 



m£ln 



n 



-x 2 H, 



>i:i 



-X3 



Here cc = (xi, x 2 , X3); p = (pi,P2,P3) and t are real variables; iro, Q±, U 2 , Q3 are real constants; 
v = ( n 2i n 3) £ z (^)i are auxiliary real functions, whose form is given below in each case, 
and the function S°(t) is given by 



S°(t) 



(X(t),P(t)) 



P 2 (t) 1 - 



(It. 



(50) 



For the expressions below to have a more compact form, we also use the following notation for 
trigonometric and hyperbolic functions 

z+(t) = -tg(t), z~(t) = th(i); z + (t) = cos(t), z~(i) = ch(t); 
z + (t) = sin(t), z~(i)=sh(i). 

The solutions of system in variations, equation (|49|). are normalized by the conditions 

{ aj (t),a fc (t)} = {a*(t),4(t)} = 0, |a j (t),4(t)} = 2i(%-^ii), fc,j = M. (51) 

The solutions of system (|49|) normalized by condition (|51|) determine the matrices Q 1 (t) and 
<5^(t), which are found as 

Qf(t) = diag-^mftf z ± [Of t], 0, 0}, Qf(t) = diag {0, rafl^, mQf}. 
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The vector jf(t) then takes the form 



1 



\z±[nft] 



,0,0 



where 7r are constants of integration. 

For the case considered we have = —Kx^, x^ 



fr 1 , and 



P ± (t) = (TRtnf ^[nft + c^],0,0) T , A ± (t) = (A ± z ± [nft + ^]/m,0 J 0) T J 

where i? 1 * 1 and i? 1 * 1 are constants of integration. Then, following (|39|). (|40|). and (|41|). we obtain 
the functions 



[x, t, h) 



\it±(t)\ ex p . 


h^ n 2' n 3 A^E J CO J ^1 J ^2 ' ^3 ' ^ (^1 0) 


y/mg ch 







(52) 



Here the multiindex v in general formula (|39j) is of the form f = (712,713) £ T? + . 

The functions (|52|) can be considered excited states for equation © with equation (|18|) being 

the linear associated Schrodinger equation. 

Note that for n 2 = 713 = 0, the expression (fB^|) gives the "vacuum" solution (|35jh 

Let (/?(x2,^3) be a function of 5^ class. Consider the Cauchy problem for equation (|4"T|) 



F0 



exp (i Im ln[(/?(x2,x 3 )]) 



=a x/m^ch (^vr (xi -xio) +Reln[93(x 2 ,x 3 )]) 



Let t^g(i) denote the linear evolution operator 



ip(x 2 ,x 3 ,t) = Ufi> o (t,s)(p(x 2 ,x 3 ) -- 
with the kernel 

00 

G$ (x 2 ,x 3 ,y 2 ,y 3 ,t,s) = ^2 

ri2=n 3 =0 



Gh ( x 2, x 3 , 2/2, 2/3, t, s)ip(y 2 ,y 3 )dy 2 dy 3 



1 



n 2 \n 3 l \2 



H, 



1)2 



Ax 2 



x H 



x exp 



AX3 ffn, 



1 



Ay 2 iT n 



Ay 3 



1 



x exp (0 2 (A^ + Ayf) + n 3 (Ax| + Ayf)) 

A Vj = Vj ~ x j0, Axj = xj - X j0 (t), j = 2, 3. 
Using the Mehler formula 



m 2 02^3 



7T 



2 h 2 ' 



00 1 

Err U ^(^"(f) 



n=0 



exp 



2xyA- (x 2 + y 2 )X 21 

H^A 2 " 



where A is an arbitrary complex parameter, |A| < 1, we obtain 

G% Q (x 2 ,x 3 ,y 2 ,y 3 ,t,s) 
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Vdh{t - s) 



cxp 



Here 



g h (t - s) = - 
f(t,s,Q,x,y) 



- (f(t, s, Q 2 , Ax 2 , Ay 2 ) + f(t, s, Q 3 , Ax 3 , Ay 3 )) 



4vr 2 /i 2 sin(0 2 (i - s)) sin(0 3 (t - s)) ' 

mn(2xy - (x 2 + y 2 ) cos(J7(t - s))) 



2sin(0(t - s)) 

Define the operator by its action on the function (/?(x 2 ,x 3 ) £ S l h as 



7T, 



where 



S±(x,i) = 5 0± (t) + P 1 ± (t)Axf + (Axf) 2 z ± [nft\ - him (ln[^(x 2 , x 3 , i)]) 



7Tn 



± 



£ In jz* [fif t] | + ft Re (ln[<^(x 2 , x 3 , *)]) • 



(53) 

(54) 
(55) 



'^'^ z±[fift] 2 

Here 5 0± is given by (JHOJ) with P(i) = P(t) ± and = X(t) ± , respectively. 

We call the operator of the form (|53|) . (|54|) . and (|55[) the "transverse" evolution operator 
for equation (|47|). To construct localized solutions for equation (|47|) . we can take, for example, 
the function 



( p{x2,x 3 ) = exp 



■^{4 + 4) 



as the initial function, i.e., at t = 0. 

Below we put m = l, g = l, K = E for = ^qo, and K = —E for = ^qq, where E is 
the unit matrix and is the "vacuum" function 1)52(1 . 

Then, for the (l + l)-dimensional case, the module of function (|52[) under the initial conditions 
P(0) = 0, x(0) = 0, vr = 1, <Ti(0) = take the form 



\^ + (x,t,, 



cos(£)| -1 ch -1 



n 



h \cos(t) 



In I cos(t) 



for 



K = k = 1 



or 



\W-(x,t,h)\ = |ch(*)j _1 ch -1 



x h 

—— In chit ) 

h \ch(t) 2 v ; 



for K = k 



-1. 



(56) 



(57) 



In Fig. [2k, the function (|56|) is shown for the (1 + 1) space-time and the potential well 
(kx 2 /2, k > 0) of the oscillator form in (|46|) with h = 1. It can be seen that the solution, 
being localized at the initial time, is focused within the time interval (0, vr/2) (in the process of 
evolution) and collapses. The GPE with an external field was used to describe the dispersion 
and diffraction of nonlinear waves [22] . The collapse phenomenon was studied in jS2] by numer- 
ical simulations. The numerical results qualitatively correlate with those obtained by analytic 
asymptotic methods. The asymptotic solution (|53jl describes the system for any T within a fi- 
nite time interval [0, T\. If the time interval is short enough, the evolution can end before the 
collapse occurs. 

The function (|56[) is -/r-periodic and its graph is given in Fig. Eb- The term (—h/2) In | cos(i)| 
in the argument of hyperbolic cosine in equation (|56|) shifts the function centroid (|56[) within 
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Figure 3. 



the time interval < t < tt/2. Throughout the evolution time, this gives rise to oscillations of 
the function maximum about the plane x = 0. 

For the potential hill of the oscillator form (kx 2 /2, k < 0) in equation (|4tij). the state of 
the system in the (1 + 1) space-time is described by ^f~. Let us choose the initial condition 
for the function the same as for Then the dynamics of the system is characterized by 
defocusing and exponential damping (see Fig- EI) ■ The term (—h/2) ln(ch(t)) in the argument of 
hyperbolic cosine in 1)5 7jl shifts the function centroid in the negative direction of x. 

Consider the (2 + l)-dimensional case. For the functions (|52|1 put P(0) = 0, x(0) = 0, ttq = 1, 
o"i(0) = 0, then we have 



ft) I 

\y-(x,t,h)\ 



cos(i)| _1 ch _1 
ch^l^ch- 1 



1 

h V^COS^) 



+ - - In I cos(t)| 



+ 



h \ch(t) 2 2 



lnch(t) 



(58) 
(59) 



The functions (|58|) and (|59() coincide at t = 0. The initial graph of the functions ()58|) and (I59j) 
is presented in Fig. Figs. \$k>,c show the graphs of function (|58j) for t = 1.2 and of the 
function ()59|) for t = 2, respectively. 

The solution considered is localized about the parabola with the vertex at the point P (t) = 0, 
X (t) = of phase space. In the course of evolution of |^ + |, the parabola branches disperse and 
at t = 7r/2 it transforms into a straight line. The inverse behavior of the parabola takes place 
for the evolution of |^ _ |: the parabola branches converge for an infinite time. The amplitude of 
the functions l^^l behaves, in the direction along the coordinate x, as in the (l + l)-dimensional 
case. 
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Figure 4. 



6 Conclusion remarks 

The asymptotics obtained, equation (JBJ), can be regarded as a necessary step in the construction 
of a global semiclassical asymptotic solution to the GPE |2j). It may be supposed that the 
functions © describe the behavior of an element of the global solution in the neighborhood 
of a ray along the normal tt (Fig. ^) to the (closed) surface about which the global solution is 
concentrated. Substantiation of these global asymptotics for finite times t G [0, T], T = const, is 
a special nontrivial mathematical problem. This problem is concerned with obtaining a priori 
estimates for the solution of a nonlinear equation, which are uniform in parameter h G]0, 1] and 
is beyond the scope of the present work. Note that, in view of the heuristic considerations given 
in [E], it seems that the difference between the exact and the constructed formal asymptotic 
solution can be found with the use of method developed in [Tal 123] . 

The technique of construction of semiclassical asymptotics developed in Section 0] provides 
a way for solving the problem of correspondence between the classical and quantum results for 
quantum systems described by nonlinear equations, namely, via finding solutions to the dynamic 
system (|26j) and (j28|) . For nonlinear systems this problem differs from the relevant problem in 
the linear quantum mechanics. For linear quantum systems, a correct formal transition from the 
quantum theory to the classical one requires imposing special limitations on the quantum states 
(their semiclassical concentration). The states not satisfying these limitations are regarded as 
"essentially quantum", and those satisfying them are considered "near to classical". Thus, the 
dynamics of the classical objects obtained is described by the classical Hamiltonian equations no 
matter the domain where the wave function is concentrated (whether it be a point, a curve, or 
a surface). The study of the dynamic system (|26j) and (|28|) is a separate mathematical subject 
for research. For example, the Hamiltonian or Poissonian formalisms as applied to this system 
are of interest. 

The construction of "transverse" evolution operator given in Section allows one not only 
to obtain semiclassical asymptotics but also to construct approximate symmetry operators of 
special type acting in the class of S l h functions under consideration. Such operators can be 
naturally referred to as semiclassical symmetry operators |33j (see also |12j). 
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